Collisionless plasma shocks in striated electron temperatures 
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The existence of low frequency waveguide modes of ion acoustic waves is demonstrated in magnetized 
plasmas for electron temperature striated along the magnetic field lines. At higher frequencies, in a band 
between the ion cyclotron and the ion plasma frequency, radiative modes develop and propagate obliquely 
to the field away from the striation. Arguments for the subsequent formation and propagation of elec- 
trostatic shock are presented and demonstrated numerically. For such plasma conditions, the dissipation 
mechanism is the "leakage" of the harmonics generated by the wave steepening. 

PACS numbers: 52.35.Tc, 52.65.-y, 52.35.Fp 



Formation of shocks as described by Burger's equation 
lH]] can be understood as a balance between the energy in- 
put by an external source (a piston moving with velocity U, 
for instance) and viscous dissipation, with kinematic vis- 
cosity coefficient i/. In one spatial dimension, this nonlin- 
ear problem can be solved exactly by a Cole-Hopf transfor- 
mation to demonstrate, for instance, that the shock thick- 
ness varies as ~ ly/U with the basic parameters of the 
problem. In principle. Burger's equation can apply for any 
continuous viscous fluid media, also plasmas. Experiments 
performed in the strongly magnetized plasma of the Ris0 
Q-machine [2] demonstrated that for moderate electron to 
ion temperature ratio T^/Ti, the strong ion Landau damp- 
ing prohibited the formation of shock. For large tempera- 
ture ratios, the ion Landau damping is reduced, and there is 
a possibility for forming steady state nonlinear shock-like 
forms, propagating at a constant speed flisll- 

In the present study we present a novel mechanism of 
an effective energy dissipation; selective radiation or "leak- 
age" of short wavelength ion sound waves. We also demon- 
strate that electrostatic shocks can form as a balance be- 
tween these losses and the standard nonlinear wave steep- 
ening as described by the nonlinear term in the "simple 
wave" equation, du/dt + udu/dz = 0, [T,*?]. Studies in 
two spatial dimensions are sufficient for illustrating the ba- 
sic ideas, and the analysis of the present paper is restricted 
to 2D. 

Magnetized plasmas are considered here for conditions 
where the electron temperature Tg varies in the direction 
perpendicular to an externally imposed homogeneous mag- 
netic field Such conditions occur often in nature for 
plasmas out of equilibrium fvl]. For the present analysis 
it is essential that the ion cyclotron frequency is smaller 
than the ion plasma frequency, i.e. ficz < ^pi- The rele- 
vant frequencies are assumed to be so low that an inertia- 
less electron component can be taken to be in local Boltz- 
mann equilibrium at all times. We assume quasi-neutrality, 
Uf, ~ rii. For a linearized fluid model we readily derive a 



basic equation in the form 



Q2 



^ = 0, 



(1) 



where is the electrostatic potential, related to the rela- 
tive density perturbations as etp/Tf, = rj = ni/ng. Since 
we are here only interested in cases where T^/Ti ^> 1, 
we took Ti = in ([T]). In case = constant, a linear 
dispersion relation is readily obtained from ([T]) by Fourier 
transforming with respect to time and space. This disper- 
sion relation contains two branches, one for uj < Q^i and 
one for 0,ci < w < Op^, the latter containing also the ion 
cyclotron waves. The wave properties of the two branches 
are very different, as illustrated best by the angle between 
the group velocity and the wave- vector [6]. For very low 
frequencies, uj <^ 0,ci, these two vectors are almost per- 
pendicular, while they are close to parallel when u ^ Alex- 
in the limit — t- 0, the dispersion relation reduces to 
(uj^ — 17^-) ((^^ — k'^Cl) = containing ion sound waves 
and the electrostatic ion cyclotron resonance. 

If we let Te = Tf,{x), with z along the magnetic field B 
and X in the transverse direction, we can still Fourier trans- 
form with respect to time and the 2; -direction. We denote 
the Fourier transformed electrostatic potential by tp. Nor- 
malizing frequencies and lengths so that 17 = w/rJci and 
^ = xQci/Cs, respectively, we readily obtain the expres- 
sion 



To 



(2) 



where we introduced a normalized propagation speed 7^ = 
(u;/fc^)^(M/To), and Tq is a reference temperature. Thus, 
7 measures the ratio between the phase velocity along the 
magnetic field and a sound speed, so that 7 = constant 
would correspond to exactly non-dispersive wave propa- 
gation. The expression ^ has the form of an eigenvalue 
equation, with 1/7^ being the eigenvalue. We present nu- 
merical solutions in Fig. [T] In the low frequency limit 
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FIG. 1. Illustrative examples of numerical solutions of ^ for 
the electron temperature profile shown in the upper panel, and 
defined by Te{^) /To = 1- ^V + V exp{-^'^ /W'^), where W = 4 
and V = 24/23. Waveguide mode solutions for Q, = 0, ^, ^, | 
are shown in the middle panel. Free modes solutions obtained for 
^ , ^ , 2 and 7 = 2 are shown in the lower panel. 
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w < rJci' shown in the middle panel, the waves are con- 
fined to the electron temperature striation (here denoted 
"waveguide modes"), corresponding to a discrete set of 
eigenfunctions ip^, with mode number m. For the Gaus- 
sian profile T^{x) considered here, the mode number m 

corresponds to the number of zero-crossing of V'm(a^) B- 
For the three dimensional problem we would have two in- 
dexes "il^km corresponding to the two directions perpendic- 
ular to B. From ip^^ we can obtain the corresponding 
eigenmodes for the S-parallel velocity u\\. The value of 
7 depends on and we find for = 0, i, i, |, that 
-il = 1.3942, 1.3805, 1.3311, and 1.2006, i.e. a relatively 
weak variation of 70 with f]. We note also that the eigen- 
functions change only little in spite of the large change in 
il. Only for close to unity, say around 0.9 or larger, do 

we see significant variations in ipoix) and 79. For shal- 
low temperature variations and narrow temperature ducts, 
we have only the lowest order mode ^^0(2;). For 7 smaller 
than the minimum value of Te{x)/To we have a contin- 
uum of eigenvalues with corresponding eigenfunctions. In 
all cases we used To = |(Te(0) + Te(|oo|). 

For uj > 0,ci the right hand side of (|2| changes sign, and 
the nature of the eigenmodes changes as well, to become 
free modes as seen in the lower panel in Fig. [T] If we let 
the electron temperature striation vanish to have a uniform 
Tg, then the free modes degenerate to two obliquely prop- 
agating plane waves. 



We consider now the low frequency limit of the branch 
of dispersion relation with oj < Qd- For m > 1, the 
waveguide modes can decay for one m-value to modes 
with other m-values. The m = mode has no decay to 
other forward propagating modes, and will be the one con- 
sidered here. For this highest phase velocity mode, with 
eigenmode i^aix), wave steepening will be the dominant 
nonlinearity llllil]. The nonlinear terms couple the var- 
ious modes to give products of x-modes. These can be 



expanded as, for instance, ipi{x)ipj 



X 



where we assume that the set ipm is complete and orthonor- 
mal [8]. For the lowest order waveguide mode we have, in 

particular, ^joo = J^ao^oix)^jix-)dx- For a large class 
of relevant electron temperature profiles we can ignore all 

higher order modes, and retain only introduce here 

^0 = ipo{x)dx. For the Gaussian variations of Tq{x) 
studied here, we will have (^0 > Cj for J ^ 1> since 
'tpQ is the only eigenfunction that is positive everywhere. 
To lowest order in the present low frequency limit we have 
the relation eip/T^ = rj ~ u\\/Cs between fluctuations in 
relative density and the velocity in the S-parallel direction. 
Our arguments concerning the mode structure therefore ap- 
ply to the velocity variations as well. 

Considering the limit of time scales much larger than the 
ion cyclotron period, we find after some algebra the result 

^ + (Con„±C.C.o)^ 



0, 



(3) 



with u\\ = u\\{z,t), where the numerical value of C,t is 
that of 7o in the limit of O — > 0. Small polarization 
drifts _L B are ignored. We introduced a constant refer- 
ence sound speed Cso- We anticipate that is here not 
much different from Qo^ since the form of ipo{x) is close to 
Tg (x) — Te ( 1 00 1 ) . The solutions of ^ have the well known 
steepening of the initial condition. The characteristic time 
for wave breaking is approximately >C/max{u||(t = 0)}, 
where C is the characteristic scale length of the initial 
perturbation along B and max{?i||(t = 0)} is the maxi- 
mum value of the initial velocity perturbation. The model 
equation ^ assumes tpoi^) ^nd 70 being used also when 
i7 > 0, but this approximation is acceptable for at least 
< < 0.75, as seen in the middle panel in Fig.[T] For 
the basic ideas outlined in the present work, this restriction 
is of little consequence. Polarization drifts become increas- 
ingly larger as 17 is increased, but these do not affect the 
dynamics parallel to B, which is covered by 

If we initialize the system with characteristic wave- 
lengths corresponding to frequencies uj <C i-S- >C ^ 
Cs/^ci, the short time evolution will be governed by (O, 
and we will have shorter and shorter scales developing as 
for the usual breaking of waves |[lll3]- This process is how- 
ever arrested when the characteristic length scales become 
of the order of the effective ion Larmor radius Cg/^d^ 
where the modes become radiating, and are no longer con- 
fined to the waveguide. We propose a phenomenological 
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FTG. 2. Examples of numerical simulations showing the varia- 
tions of the normalized potential $ = etp/Ti as a function of po- 
sition at a fixed time. The electron temperature enhancement is 
localized as a Gaussian in the x-direction. We have uj = Q,pin/5 



0.05 andQ,ci/0. 



in both cases while Qci/^pi 

and right panels respectively. We have here Te/Ti — 50. 
clarity, only a part of the simulation domain is shown. 
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expression for the process, best written in Fourier space in 
a frame moving with C<;o, as 

du 

" " ■ n 



'II , • CoA;||^ ^ 



^ 1 , (4) 



with U|| = 'U||(/cj[,t). The symbol ® denotes the convolu- 
tion product and % is Heaviside's step function. T char- 
acterizes the time it takes for the energy of the uj > Qci- 
waves to be lost from the waveguide. We have T = T{k\\), 
but it will depend also on parameters such as the waveguide 
width as well as the other plasma parameters. We expect 
that increasing width gives increasing T, i.e. decreasing 
shock thickness A. Within the present model the waveform 
will steepen uninhibited until the shock thickness becomes 
of the order of Cg/Vt^i, at which time the harmonic fre- 
quencies will exceed Vt^i to become radiative and the high 
frequency wave energy is lost from the waveguide. Conse- 
quently, the S-parallel scale of the shock is controlled by 
a quantity referring to the S -perpendicular dynamics. 

We study the nonlinear propagation of low frequency 
waves in an electron temperature striated magnetized 
plasma numerically by using a 2 i -dimensional particle in 
cell code described elsewhere |@]. The code assumes ex- 
plicitly the electrons to be locally Boltzmann distributed 
and the resulting nonlinear Poisson equation is solved by 
iteration. We use a Gaussian variation for T^{x) so that 
TedooD/T, = 1, while TM/T^ > 1, where T, = 
constant. The width of the electron striation is here 
W = 28Ai5, = 3.9C,oMi- Several values of re(0)/r, 
were investigated. 
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Fit result: 5n/n =0.24±0.002, z =824±1 A=70±l 
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FIG. 3. Numerical simulation showing normalized potential 
$ = eil}/Ti and relative density ni/riQ during the formation and 
propagation of a shock under the conditions mentioned before. 
We have Te/Ti — 25, ^^.i/^pi = \ and Sni/riQ « 0.24, where 
Srii is the actual detected density perturbation at the time where 
the shock is fully formed. The background density uq is normal- 
ized to unity. Lower frame shows sample of shock fitting for the 
ion density. 



Results illustrating the waveguide and the free modes are 
shown in Fig. |2l The properties are clearly different and 
consistent with the interpretation given before. The wave- 
guide modes are confined to the electron striation, while 
the high frequency free modes are dispersing or "leaking", 
consistent also with laboratory experimental results [S]. It 
is important to emphasize that this apparent damping will 
be found also in a fluid model. The observed effective 
damping is caused by wave energy dispersing in space, 
and dissipated by linear ion Landau damping outside the 
electron temperature striations where Te(joo|)/Tj = 1 in 
all cases considered. Each step in the process is formally 
time -reversible. In order to emphasize the physical effects 
we discuss here, we consider only high temperature ra- 
tios, Te(0)/Ti > 25, in order to reduce the effects of lin- 
ear as well as nonlinear ion Landau damping. Such high 
temperature ratios (even as large as T^/Ti = 100) can 
be obtained in discharge plasmas under laboratory condi- 
tions ||9|]. For nonlinear waves described by a Korteweg-de 
Vries equation, a shock-like structure is followed by Airy- 
type ripples, originating from the dispersion term in the 
KdV-equation. These ripples are absent in our results, see 
Fig. |3] Likewise, at the high temperature ratio used here 
(with Cso ^ Uth.i, the ion thermal velocity), we find no 
ions being reflected by the shock. A backward propagating 
rarefaction wave is of no concern here. 

In Fig.[3]we show an example of the formation and prop- 
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FIG. 4. Shock velocity U in units of Cso (upper panel) and sat- 
urated shock thickness A in units of Cso/^ci (lower panel) as 
function of the fully formed shock relative height Sni/riQ, for dif- 
ferent combinations of the parameters Te /Ti = 25 and 75, and 

^ci/^pi ~ 5 and 2' 



sponds to approximately 50 % reduction in A. We demon- 
strated also that the shock thickness A scales inversely pro- 
portional to the width of the striation. We verified that 
changes in the initial S-parallel scale length of the initial 
density profile do not change the saturated shock thickness 
A for any Srii/riQ. We studied also a weak magnetic field 
limit, with Qci/^pi = j- Here, finite ion Larmor radius ef- 
fects (ignored in deriving (dJ) begin to be important, since 
the striation is now only ~ 10 ion Larmor radii wide. For 
this limit we find that A is almost constant. Even weaker 
magnetic fields will require a fully kinetic theoretical anal- 
ysis to account for the effects of collisionless ion viscosity 

We here reported arguments for the formation of shocks 
in electron striations in magnetized plasmas when Q^i < 
Qpi and Tg ^ T^. By PIC-simulations we demonstrated 
the formation and propagation of electrostatic shocks un- 
der these conditions for a wide combination of parameters. 
The decreasing shock thickness for increasing amplitudes 
is found also for classical shocks as described by Burger's 
equation, but in our case the dissipation mechanism is leak- 
age from the electron temperature striation of the short 
scale lengths generated by the nonlinear wave steepening. 
The ideas presented here can have wider applications. 



agation of a shock. The initial density has an error function 
type spatial variation in the S-field aligned direction. Ions 
are continuously injected at the boundary at z = to main- 
tain the energy input. We have T^/Ti = 25, ^Ici/^pi = \ 
and bnijriQ ~ 0.21, where brii refers to the actual de- 
tected density perturbation at the time where the shock 
is fully formed, and not to the initial imposed perturba- 
tion. The standard deviation of the parameters is estimated 
from fitting the shock profile with a nonlinear Levenberg- 
Marquardt method. 

In Fig. |4] we show the shock thickness A and the nor- 
malized shock velocity U for various combinations of the 
parameters T^/Ti, ^^i/^pi and dni/riQ. The velocity U 
is obtained by following the shock position as func- 
tion of time. For small or moderate values of Sni/riQ we 
find a close to linear relationship between 6ni and U. As 
6ni — 7- (shown with dashed line fit) we have U approach- 
ing the appropriate value for 70, apart from small correc- 
tions due to finite Ti ignored in dJll. For Qd/^pi = 5 we 
find that an approximate five-fold increase in 6ni /uq corre- 
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